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1. Introduction 

The deflection of photons from distant astronomical sources by intervening massive 
objects is rapidly becoming an important tool for experimental astrophysics and 
cosmology. At the same time a variety of optical, atomic and molecular interferometers 
in use or under development is approaching the threshold of measurability for inertial 
and gravitational effects. 

Photons in gravitational fields are described by the generalized Maxwell equations 



where V Q indicates covariant differentiation. We use units h~ = c = 1. 

In interferometry, where the phase shifts are proportional to the mass of the particle 
employed, one may wish to consider, beside photons, also massive, charge-less, spin-1 
particles that satisfy the generalized Proca equation 



in the hope that among the atoms and molecules to be used some indeed obey (J3]). 

We have shown in previous work p]-[6], that equations ([I])-© can be solved exactly 
to first order in the metric deviation 7^ = — r]^ u , where l^^l *C 1 and 77^ is the 
Minkowski metric of signature -2. In all these solutions, the background gravitational 
field appears in a phase operator that acts on the wave function of the relative field-free 
equation. It is precisely from the corresponding gravity-induced phases that the optics 
of the particles considered follows. 

The purpose of this paper is to show that the use of the solutions found presents 
advantages, in principle, over other current approximations. This is particularly true of 
gravitational lensing where large use is made of the thin lens approximation [7], [8] in 
which the lens is characterized by a surface mass density and the waves are scattered 
only at the thin lens plane. The phases can in fact be calculated by means of path 
integrals that can be integrated easily and exactly for most known metrics. Gravity 
acts along the complete length of the space-time paths and not only in a narrow region 
at the lens plane [9]. 

The plan of the paper is as follows. In Section II we give, for completeness, the 
solutions of (JT])-(J3|) and show, in Section III, that they are gauge invariant. Geometrical 
and wave optics are discussed in Section IV. We summarize the results in Section V. 

2. Solution of the spin-1 wave equation 

To first order in 7^, ([I])-© become 



V a V a A, = 0, 
or, more generally, by the Maxwell-de Rahm equations 
V a V a A lt -R fia A a = 0, 



(1) 



(2) 



V a V a A^ + m% = , 



(3) 



V,V^ ~ (7/ 



- 2T 9tia/ A^ - \l a ^A° = , 



(4) 
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V.VM, - R, a A a ~ (77™ - 7 ™)4^ + 21^,4^ = , (5) 

(V.V+m 2 )^ ~ {^-Y a )A^2Y a ^ v A^~ la ^+m 2 A^ = , (6) 

where ordinary differentiation of a quantity $ is equivalently indicated by $ )Q or <9 a <I> 
and r CTj ^ = 1/2(7^^ + 7 

— ltiv,a) <wg the Christoffel symbols. In what follows we 
also use the notations K a = x a — z a and B{ a \...m = B a \.„p — Bp\ ^ a . In deriving (jl])-©, 
we have used the Lanczos-De Donder gauge condition 

7 " - -7 CT = . (7) 

lav, 2 ^ ' 

The field A^x) also satisfies the condition 

V^" = . (8) 

Though (jHJ) is strictly required only when m = 0, it is convenient to impose it also in the 
case of a massive spin-1 particle. Equations (jlj) and (jBJ) can be handled simultaneously. 
As shown in Appendix A, their solution is 

1 f x 

A^x) ~ e'^a^x) « (1 - = ^(^O + 9 ' ^^W^/^) ( 9 ) 



p 



p Jp 



where a M satisfies the equation d v d u a^ = in the case of (jl]), and {d u d u +m 2 )a^ = when 
flSJ) is a solution of (jHJ). The solutions a M (x) used in this paper (Section IV) are either 
plane waves or spherical waves multiplied by e~ lk ° x ° . For this reason, it is convenient to 
group % with £. 

Similarly, it can be verified by substitution that the Maxwell-de Rahm equation (jSJ) 
has the solution [I] 



A^x) ~ e « (1 - irfia^x) =a fl + j dz ^(z)d x a a (x) - (10) 



where again d a d a a fl = and 

- irja^x) = -i£a p (x) - - J dz x [y a ^ x (z)a a (x) - 2'j a ^(z)d x }a a (x) . (11) 

The last two terms are typical of the Maxwell-de Rahm solution and have been tagged 
with the parameter e (e = 1) for bookkeeping purposes. Solutions (jHJ) and (jTDJ) must 
also satisfy (jSJ). On using V a A a ~ rf-® A a $ — 7 a/3 a Qj/ 3 and (1A.3J) . we obtain 



V a A Q = f(x) (12) 

1 



dx x 

-^-7 af ,(x)d^dxa x (x) 
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By performing the gauge transformation A a — > A a -\- A )(T , we find V a A a — > V a A a + 
V^Vo-A + / = 0, from which we obtain V cr A a = provided 

V CT V (7 A + / = 0. (13) 

Finally, the addition of a term rn^A^ to the l.h.s. of ([2]) invalidates the solution (jTOl) 
that can not, therefore, be extended to include de Rahm's term, except when 2em 2 ^ iiv a v 
can be considered a second order quantity. It is known, on the other hand, that the 
term R^ a A a can be dropped in lensing when the wavelength A of A a is much smaller 
than the typical radius of curvature of the gravitational background, or when, in the 
weak field approximation, R^ a = — l/27„ a ^ = 0. 

3. Gauge invar iance 

The first two integrals in (Q and ([TO]) represent by themselves a solution of the Klein- 
Gordon equation (V^V 1 ' + m 2 )0 = 0. The additional terms must therefore be related 
to spin. In fact (Q can be written as 

A^ = cif_i — — J dz x 1[am K a d p a^ + -J dz x laX d a a^ + (14) 

+~J dz\ am S af3 a^+ l -j^ dz^xT^a^x) , 
where the tensors 

(s a "U = - <W > ( TQ/ V = + (is) 

have been written in matrix form as S a/3 and T af3 and S al3 a fJ- = (S a ^) ^a 1 ' , where a v 
is a column matrix. The matrices Si = 2eijkS^ k satisfying the commutation relation 
[Si, Sj] = ieijkSk, can be recognized as rotation matrices. These, and the matrix a u , are 
represented by 



Sr = - 



s 3 = - 
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We give the matrices S a ^ in Appendix B and show the equivalence of matrix and tensor 
expressions. In the same Appendix we also show that by applying Stokes theorem to 
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the line integrals on the r.h.s. of (I14p extended to a closed space-time path I, we obtain, 
to first order 1101. 



A " = ( 1 ~ \ 7 S d7 " SR °W °M > ( 16 ) 

where J a/3 = L a/3 + S a ^ is the total angular momentum of the spin-1 particle, 
Rtivap = l/2(lfni3,va + lva,nP ~ lfj,a,u/3 - Ju/3,»a) is the linearized Riemann tensor, and 
E is the surface bound by T. Equation ffT6l) is unambiguous because, to first order, the 
values of 7^ on S can be obtained from those on T independently of the path followed 
[TO] . The term of (TH|) that contains S a p gives rise to the Skrotskii effect [llj,[12j. It 
follows from (TIT)]) that ([2]) is covariant and also invariant under the gauge transformations 
7iw ~~ > 1/j.u — i^,v — £,v,n induced by the coordinate transformations x M — > x^ + £ M which 
are still allowed in the weak field approximation. The quantities also small of first 
order, must satisfy the equation £ ° = because of the Lanczos-De Donder condition. 
It also follows from (fl6|) that the term containing T af3 in (I14p does not contribute to 
integrations over closed paths, behaves as a gauge term and may therefore be dropped. 

The difference between ( !T0|) and (Q is represented by the two terms labelled by e 
in ( TlTl) . Of these, the first one vanishes when integrated over a closed path. To evaluate 
the second term we apply the transformations 7^ — > 7 M „ — — ^ jjU , integrate over a 
closed path and choose, for simplicity, £ Q = (, a oe~ dx and a a = a a0 e~ lkx , where l a l a = 0, 
k a k a = 0. We obtain 

-ie j> dz x k x {£ a ^{z) + £^ a {z))a a {x) = -^l[^ a }oa a (x) j> dz x l x + 



+ie 



1,, (b dz (d x a a (x)^ a0 ) i> dz (d x a a (x)l a 







because both integrands are of the type du and return to the initial value at the final 
point. This shows that the addition of the de Rahm term to the solution does not affect 
its invariance under the gauge transformations of 7^. 

We must now consider the effect of the electromagnetic gauge transformations 
o>fi — > + A iAt which have already been used to obtain ( Tl3l) . Equation ([9]) can 
be explicitly re-written in term of the gauge invariant electromagnetic field tensors 
F^v = V M A, - VuAf, = A M , and = a M as follows 



F, v — fnu 2 i dz\ aXA {z)k a d f, u {x) + l_ dz\ x m (z)U(x) 



p J p 



+1 I dz x laX (z)d a f^{x)+^(x)J rS[ r) 



p 



The additional terms that the Maxwell-de Rham solution ffl0|) has relative to (Q can 
be re-written in the form 

' dz x [y^ x ( Z y(x) - 2rf(z)f Xa (x) - 2 1 «(z)a x , a (x)] . (17) 
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The second term of (IT7|) is already invariant because it contains f\ a . The first term 
vanishes when the integration is over closed space-time paths. The last term transforms 
according to 

/X PX 
dz x ^(z)ax,a(x) e J dz x ^(z)(a XiCt + A M ) , (18) 

can be re-written in the form 

and vanishes when the integration is over closed space-time paths. Equation (TTUT) is 
therefore invariant under electromagnetic gauge transformations. 




4. Optics 

4-1. Lensing 

The correct general relativistic deviation of light rays in the gravitational field of a 
star follows immediately from the first two integrals in (Q and ({TO]) . We choose, 
e.g., a gravitational background represented by the Lense-Thirring metric, 700 = 20, 
7ij = 2(f>5ij, (f> = —GM/r, and 7 j = hi — 2GJijX^/r 3 , where x J = (x,y,z), 
r = \J 1 x 2 + y 2 + z 2 and Jij is related to the angular momentum of the gravitational 
source. In particular, if the source rotates with angular velocity u = (0,0, a;), then 
h x = 4GMR 2 ujy/5r 3 , h 2 = -4GMR 2 ujx/5r 3 , . 

Without loss of generality, we assume that photons propagate along the z-direction, 
hence k a ~ (k,0,0,k). Moreover, photons propagate along null geodesies, from which 
it follows that ds 2 = 0, or dt = dz. Using plane waves for a M (x) = a°^exp{—ik a x a ) , the 
solution has the form = e~* x a°, where 

X = Kx a - ~ jf dz x 1[aK(3] {z)K^ + X - jf dz x laX {z)k a . (19) 
We can define the photon momentum as 

k a = P^- (20) 

dx a v ' 

It is easy to show that x satisfies the eikonal equation g al3 x,aX,/3 — 0- We can also prove, 
by direct substitution, that the inclusion of de Rahm's terms (see fTTTT) ) in f|T9|) does not 
invalidate the eikonal equation. Geometrical optics is not therefore affected by the term 
R^A a in ©. 

For the Lense-Thirring metric, \ is given by 

k r Q 

X - ~ 2 / {(* - x')<t>,z'dx' + {y- y')<t>, z <dy' - (21) 



p 



- 2[(x - x')<f> >x i + (y- y')4> tV >]dz'} + k j dz 
k 



Q 

{[{x - a/)(hi tZ > - 2h 3jX f) + {y- y'){h 2 , z > - 2h 3iy >)] dz'- 
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- [(ar - x')hi jX > + (y- y')h hy ,] dx' - [(x - x')h 2 , x > + (y - y')h %y >] dy'} 
k r Q 

+ - / [2h 3 dz' + hidx' + h 2 dy'} , (22) 



2 . , 

where P is the point at which the photons are generated, and Q is a generic point along 
their trajectory. The components of the photon momentum are therefore 

~ kl = 2k J P {-\f z dx -Yt dy+ d -^^ dz 

- \{HQ) ~ hi(P)) , (23) 

+ ~(h 2 (Q)-h 2 (P)), (24) 

~k 3 = k(l + <f) + h 3 ) . (25) 

We also have k = kj_ + k 3 e 3 , kj_ = k\ e x + k\ e 2 , where kj_ is the component of the 
momentum orthogonal to the direction of propagation of the photons. 

Since only phase differences are physical, it is convenient to choose the space- 
time path by placing the photon source at distances that are very large relative to 
the dimensions of the lens, and the generic point is located along the z direction. We 
therefore replace Q with z, where z ^> x,y. Using the expression for h\ 2 we find that 
their contribution is negligible and f[23]) - (T25|) simplify to 

k = 2k r gjj+M^ (26) 

J-oo dx 

~ k2 = 2k J*JJt±M dz , (27) 

~k 3 = k(l + (p + h 3 ). (28) 

From (!26l) - (f28l) we can determine the deflection angle 9. Let us analyze first the case of 
non-rotating lenses, i.e. h 3 = 0. We get 
~ , 2GM ( z\ 

k 1 ~h n p-x(l + -), (29) 
~ , 2GM ( z\ 

k 2 ~h—— y 1 + -) , (30) 



R 2 ° \ r, 

~k 3 = k(l + <p + h 3 ), (31) 



where R = ■sjx 1 + y 2 . By defining the deflection angle as tan^ = k±/k 3 , it follows that 

„ „ 2GM ( z\ 
tan 9 ~ 6 — 1 + - . (32) 



R \ r, 

In the limit z — > oo we obtain the usual Einstein result 9m ~ 4GM/R. The remaining 
term of (TH|) contains S a p and therefore is a spin effect. It is usually referred to as the 
Skrotskii effect and represents a rotation of the plane of polarization of the incoming 
photon beam. 
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A general expression for the index of refraction n can also be derived from (1191) . (1201) 

and n = k/k = x,a/Xfi- 

4-2. Wave effects in gravitational lensing 

We consider the propagation of light waves in a background metric represented by 
7oo = 2U(p), 7ij = 2U(p)Sij, where U(p) = —GM/p is the gravitational potential 
of the lens. Wave optics effects can be calculated by using the type of double slit 
arrangement indicated in Fig.l. We neglect all spin effects. This limits the calculation of 
the phase difference to the first two terms in ([9]) and ( jTUl) . We also assume for simplicity 
that k 1 = 0, so that propagation is entirely in the (x 2 , x 3 )-plane and the set-up is 
planar. We use a solution of the equation <9 7 <9 7 a a = of the form a M = a°e~* fc<TX ' CT /r = 
a o e -ik x° e -ikix % J r _ a Q e -ifco*V fc 7r, valid everywhere outside the light source (r = 0), 
which is the region of physical interest in lensing. We neglect the contribution from the 
de Rahm term because, for the metric used, R ail A a = —l/2j a "A a ' = —V 2 {GM / p)A jl 
which vanishes for p ^ 0. 

The corresponding wave amplitude <fi is therefore 

<f>(x) = e -^— 1 2 + J [cf A 00 n° + dz 2 l22 Ii 2 + ^ 3 733 n 3 ] - (33) 
[dz° (7oo, 2 ^ [0 n 2 l + 7 00i3 ^[0n 3 l) + dz 2 l22 , 3 K^ + dz^ 2 K^] 



where II = — ik, IP = — ik l — x l /r, and we have taken into account the fact 
that 7n plays no role in the planar arrangement chosen. The change in phase 
must now be calculated along the different paths SP+PO and SL+LO according to 
A(f) = (0slo — 4>spo) I (e~* fc<Tl ' CT /r) and taking into account the values of II* in the various 
intervals. The total change in phase is A0 = A(psi + A(j>LO ~ Acfisp — A0po- The 
details of the calculation of A0 along the path segments are given in Appendix C. All 
integrations in ( 1C.2I) . (1C.4|) . (1C.6I) and (1C.8I) can be performed exactly. All results can 
be expressed in terms of physical variables r s , r , b + , b~, and s and lensing variables 
D s , D ds , D d , 6 + , 6~, and /3. We obtain 

A(p = y | In (^-\J D 2 dS + (s + b~f + b~ cos 7 + r^j - In (b~ (1 + cos 7)) 

+ In (b + (l — cos<^ + )) — In [rs — r^ — b + cos</? + ) 

+ In ^b~ + r o cos0~ — \Jb~ 2 + r 2 ^j — In (r (l + cos^ - )) 

+ In (r (1 + cos 0+) ) - In (V + r cos 6+ - y/b+ 2 + r$J J , (34) 

where r| = b +2 + r 2 L + 2b + r L cos^ + , r\ = D 2 dS + (s — b + ) 2 , ^ + + a + + a~ +^ — 6 + — 6~ = n 
and y = 2GMk. Equation fl34l agrees with that obtained in [9 J for plane waves and can 
be applied to the optics of any particle when spin effects are neglected. 
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Figure 1. Geometry of a two-image gravitational lens or, equivalently, of a double 
slit interference experiment. The solid lines represent the particle paths between the 
particle source at S and the observer at O. M is the spherically symmetric gravitational 
lens. S, M, O and the particle paths lie in the same plane. The physical variables are 
r s> r o> b^, s, while the lensing variables are indicated by D^g, D^o, Dg, (5. The 
deflection angles are indicated by a ± . 



The contributions % o f th e real parts of IT 2 and II 3 to (jCj) , (jClij) . (IClBj) and flCl8l) 
are listed below by path segment 

isL = S b+ + 3 erf y + , n r, - r, - 6+ cos f _ 
ri [ cos ip + b + (1 — cos(^ + J 

2r^b + sin 2 99" 
cos<^+ 



CM J 2-3cos 2 #+ 6+ - ^rfTF 2 + &+cos0+ 
R{ cos#+ r (l + cosfc> + ) 



2r v^ + & +2 sin 2 6 + 
+ b+ cos 9+ 



GM f 2-3cos 2 7l r s~ yD 2 dS + {s + b f + b cos 7 

XSP = -&-{b- -In v - — — r + (37) 

Hr cos 7 b (I + COS7) 



2b ~\/ D ls + (s + 6-) 2 sin 2 7 
rs cos 7 



GM J 2-3cos 2 fl- Zr - v 7 ^ + b~ 2 + r o cos0- 

Xpo — 1 r ° Z ln /1 1 7J^ + 3 ^ 

Ri cost'" r (i + cos6>- x 



2r a/?"q + b~ 2 sin^ 
6~ cos6>~ 
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Figure 2. The plot of F/{AGM/R X ) versus 9 has a singularity at 9 = 0.57483 where 
the argument of the logarithm in (|4"2| vanishes. For all other values of 9 and for 
reasonable values of M and R\, the amplification factor F « 1. 



Contrary to A0, x = + x LO + xsp + Xpo is completely independent of y. 

The meaning of fl35|) - fl38|) becomes evident by recalling that the total wave 
amplitude has the form 



:) = exp [— ikox — i£ (x) + x (x)] 



and that the ratio 



Wo 



3 X+X* 



(39) 



(40) 



is known as the amplification factor. As an example, let us consider the simpler case of 
particles coming from x 3 = — oo ((p + = 7 = 7r/2) with b + = b~ = s = b , r = D^s- Then 
Fig.l gives 6 + = 9~ = 9 , and R\ = R2 ■ In this case xsl = —Xsp,Xlo = Xpo and 
A(f> = 0. A simple calculation shows that the probability density of finding a particle at 
O is 



where 



X 



AGM 



x+x* 



cos 



-^e 



x+x* 



. 2-3 cos 2 9 , sin 6» - 1 + sin 9 cos 6> 
sin H : In 



(41) 



(42) 



i?i [ ' 2 "' cos^(l + cos^) 
The quantity Fj (AGM/Ri) as a function of 9 is represented in Fig. 2. It has a singularity 
at the value tan# = 0.647799 for which the argument of the logarithm in (1421) vanishes. 
We find Fwl for all other values of 9 and for all reasonable values of the parameters 
involved. 



5. Conclusions 



Covariant wave equations for massless and massive spin-1 particles can be solved exactly 
to first order in 7^. The solutions are covariant and invariant with respect to the gauge 
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transformations of and 7^ and are known when a solution of the free wave equation 
is known. 

The external gravitational field only appears in the phase of the wave function. 
We have shown elsewhere [2], [9] that the spin-gravity coupling and Mashhoon's helicity- 
rotation interaction [13], [2] follow from —i^a^ and —ir/a^. According to equations 
ffTUj) . ffTTl) and fTH|) . the spin term S a p finds its origin in the anti-symmetric part of the 
space-time connection. In the case of fermions, S a p is accounted for by the spinorial 
connection [6]. The term of that contains S a p gives rise to the Skrotskii effect. 

If the solution of the free wave equation is a plane wave, then the change in phase 
is entirely given by A<p. It is from the gravity- induced phase that we have derived 
the optics of the particles. The phase can be calculated by means of path integrations 
that can be frequently performed in an exact way, when the metric of the gravitational 
background is assigned. 

From the phases we have derived the geometrical optics of the particles and verified 
that their deflection is that predicted by general relativity. In addition, the background 
gravitational field acts as a medium whose index of refraction can be calculated for any 
metric from ( 1T91) . ( 1201) and n = k/kg. 

Wave optics can also be extracted from the gravity-induced phases. In lensing 
this is normally accomplished by means of a thin lens model whereby the gravitational 
potential vanishes everywhere except on the lens plane. In the present approach the 
gravitational field acts all along the particle's trajectory from source to observer and no 
thin lens approximation is required. 

The exact expression for the phase change A0, gives rise to interference and 
diffraction phenomena. In gravitational lensing, wave effects for a point source depend 
on the parameter y and different values require, in general, different approximations to 
the solution of the wave equation. In particular, diffraction effects are expected to be 
considerable when y ~ 1. In the present approach and in the configuration of Fig. [IJ 
( 134)) applies regardless of the value of y when A0 is real. 

The extension of our results to include spherical wave solutions of Helmholtz 
equation yields an amplification factor fHUj) whose value can be calculated exactly. We 
find that essentially F ~ 1 for all reasonable values of the parameters involved. 

Because spin has been neglected in the example given, the same results can also be 
applied to the case of gravitational lensing of gravitational waves [T5]. 
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Appendix A. 

In order to prove that ([9]) is a solution of ([6]) and therefore also of (00) when m = 0, we 
make use of the formula 

/x r-x 
dz x G x (z,x) = G T (x,x) + J dz x d T G x (z,x) . (A.l) 

By writing in the form 

A^x) ~ (1 - z£(:r)K0r) (A.2) 
and differentiating —i^a^ once according to ( 1A.1I) . we obtain 

dA-ifrJ = -\ P dz x 1[am (z)[S^a,(x) + K^ a d T d^a p {x)} (A.3) 



+^lar(x)d a a ll (x) + - I dz 'y aX (z)d T d a a ll (x) 



4 . , 

' p 

ur Al 

ip 

Differentiating (1A.3I) once more and contracting the indices of differentiation, we find 
d^d-ri-i^) = -h^dPa, - \ [ X dz x 1[rX ^d?a, (A.4) 



+T a ^ T (x)a a (x) + / dz T a ^ x (z)d T a a (x) . 



/•a: 

/ dz X 'y<,.^ m 



\ £ dz x 1[aXA K a d T d T d^a, + i 7 ; T «9 Q S + 7ar^ r S 

-j PX PX 

+ 2j P dzX ^xd T d T d a a^ + d T {T a ^ T )a° + 2T a ^ T d T a° + dz x Y a ^ lX d T d T a° 

We assume that a^{x) is well behaved and that successive differentiations applied to it 
commute. Then the second term on the r.h.s. of ( 1A.4I) vanishes because two symmetric 
indices are contracted on two anti-symmetric indices. The third, sixth and last terms 
may be written as —m 2 (—iC,)a II . We find 

ardri-itaj = -m 2 H£K - \{i a a , p - 2^)^ + 7*^% - (A.s) 

The second term on the r.h.s. of ( 1A.5I) may be dropped on account of (jZj). We now 
substitute the solution flA.2[) into the r.h.s. of ([HD and keep only first order terms. We 
get 

V.VM, + m 2 ~ {rf a - ^ a )d a d a a^ + d a d a (-iZ)a^ (A.6) 

and, on using (IA.5I) . 

V.VM, + m 2 A, ~ -i[(7^ - 2 T y ^ S + ( 7 / - 7ct ./K] = (A.7) 



2 

on account of (J7J. 
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The matrices S a/3 can be easily obtained from their definition ( 1T51) where fi and v label, 
respectively, rows and columns. We find, for instance, 

/oioo\ / o o 1 o \ 



S 



01 



V 



-10 





S 



02 



/ 



V 




-10 




Let us define the quantity A a/3 = dz x ^^ aX ^-\. On using (TO)]) , we obtain 



af3\pv 



(B.l) 



--A a/3 (77^77^ - ^V") a u = -i (A^a + A"V + A» 2 a 2 + A^a 3 ) , 

which is a four- vector and can also be understood as a single-column matrix whose rows 
are obtained by setting \i = 0, 1, 2, 3. On the other hand, using the matrices S af3 , we 
also find 



(A aP S^)(an = 

= 2 [A 01 S 01 + A 02 S 02 + A 03 S 03 + A 12 S 12 + A 13 S 13 + A 23 S 23 } 
( A Gl a l + A 02 a 2 + A 03 a 3 \ 



(B.2) 



-A 01 a° + A 12 a 2 + A 13 a 3 
-A 02 a° - A l2 a l + A 23 a? 
-A 03 a° - A^a 1 - A 23 a 2 



which coincides with the final matrix in (IB. 1[) . Therefore a u is understood as a four- 
vector in (IB.lj) . while in flB.2j) (a M ) is a single-column matrix. This proves the equivalence 
of the matrix and tensor expressions used in ( TH1) and ( |T5l) . We now apply Stokes theorem 
to 



aJx) + (t> dz G\(z,x)aJx) 



(B.3) 



where T is a closed path in Minkowski space. Following textbook procedures [10J, we 
obtain 



where 



Ap, — U.;, 



a > + 2L dT °°{£t G °-£; 



a 



I 1 

-lp[aa,(3}{z)K a d p a^(x) + -^ aa (z)d a a fl (x) + 



(B.4) 



(B.5) 



+~7 M (2)£ Q V(z) + - A l a pAz)T af3 a,{x) 
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and £ is a surface bound by T. We then find 

Grjdft = [--l[ a „,i3]8K a dP + ^[Sa,l3]d P + 



and again 



rfr CT5 (G ff) * - G 5)Cr ) ^ = - [- (7 QCTj/3 5 + 7p5,aa ~ l a 6,/3v ~ Ij3*,a5 ) (x a ~ Z a ) 8 P (B . 6) 
+ -(7ckt,/35 + 7/35,cwj — 7a5,/3<J — lf3a,as) S^a^dr" 5 = 



dr^ l -[-R ap5a {x a - z a )dP + l -R aP 8aS^\ 



2 L , 2 

By introducing the angular momentum generators of the Lorentz group 

L a(3 = (x a - z a ) {id 13 ) - {x? - z p ) (id ) (B.7) 
into (1B.6P and using (1B.4j) . we find 

A^-a^-*- f dT° s R aPaS {L<*P + S aP )a„ , (B.8) 



which corresponds to ( 1161) . 
Appendix C. 

In order to calculate A0, it is convenient to transform all space integrations into 
integrations over z°. Along SL we have 

U = ~fou/2 > ^l(z°) = (r L - z ) 2 + b+ 2 + 2{r L - z°)b+ cos ^ , (C.l) 
k 2 = k cos ip + , k 3 = k sin ip + , 

t-t2 t + coso9 + , ., . , smip + n 

II = — ik cos <p H , 11 = — ik sin H , at z = t\l , 

where sin </> + = Das- We find 

A05L ^ 



GM 



\ L dz°qsL(z y 1/2 [n° + cos^+II 2 + sin^+II 3 ] + (C.2) 
Jo 



+ 2 dz°q SL (z°)- 3 / 2 (-z + r L + b" 
Jo 

r „ 9 sin 2 oj + „ocos 2 o? + „ ni 

x -n 2 — ^— - n 3 — v + n° 

cosf + sin<£>+ 



cos Lp + )(r L — z°) x 



where II 2 and II 3 must take the values specified in (1C.1I) at the upper integration limit. 
Analogously, for LO we get 

GM 



U = n^^ 2 = ksm9 + ,fc 3 = kcos6 + ,Rx= Jr 2 + b+ 2 , (C.3) 

qLo{z u ) 1 v 



Optics of spin-1 particles from gravity-induced phases 



15 



q LO (z°) = (yjr 2 + b+* - z° + r L ) 2 + r, 



-2r (Jr% + b+ 2 - z° + r L ) cosfl H 



T-r? ■, • „4. sin6> + - ., „, cos6 l+ n 
n 2 = -ik sin 6» + H — , IT 3 = -ik cos 9 + H — , at z° = r L + R x 



Ri 



Rx 



and the change in phase is 

U <r L 



GM 



j L+ 1 dz°q~ L l J 2 [Ii + sinfl+IT 2 + cosfl+n 3 ] + 

J rr. 



(C.4) 



+ 2 / dz°q LO (z )- 3/2 (z - R 1 + r cos9+) 

Jr L 



x (r L + R 1 - z°) 



n 



,cos 2 fl+ 
sin 9 + 



osin 2 ^ 4 " « 
cos 6 ,+ 



Again, II 2 and II 3 take the values calculated in (1C.3I) at the upper integration limit. 
For SP we find 

GM 



U 



qsp^) = b- 2 + (R- z {> Y - 2(R - z»)b~ cos 7 , (C.5) 



qsp(z°) 1/2 
k 2 = k cos 7, fc 3 = /csin7, 

,R= ^D dS + (s + b 



R 2 + b -2_ r 2 



cos 7 



2b~R 

It 2 = — ik cos 7 + 



2 ■•/•-..■- i C ^H ; n 3 = -life sin 7 + , at z° = R 
R 



SU17 

IT' 



and the corresponding change in phase is given by 

A05P 



GM 



= / ^Vp(^°) _1/2 [n + cos7n 2 + sin7n 3 ] 

Jo 

+ 2 / dz°qsp(z )-^ 2 (-z° + R-b~ cos 7 )( J R 



(C.6) 



Z°) X 



X 



„ 9 sin 2 7 ~ cos 2 7 ^ n 

-n 2 - n 3 — — - + n° 



cos 7 



sin 7 



Finally, for PO we obtain 

GM 



U 



k 2 = — ksinO , k 3 = k cos 9 , 



qPoW 2 

q P o(z°) ee r 2 + (i2a + R - z ) 2 - 2r (R 2 + R - z°) cos9~ 
R 2 = sjr 2 + b- 2 , 



(C.7) 



IT = ik sin 9 — 



Sm6> „ q , „ COS# n „ „ 

— — , II 3 = -tk cos 9 H — , at z" = R + R 2 

Ri Ri 
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and the corresponding change in phase is 



>PO 



GM 



= / d Apo 2 [n° - sin rn 2 + cos r n 3 ] 

JR 

rR+R 2 

+ 2 / dz°q P o(z )~ 3/2 (z -R 2 -R + r cos 6~) x 

JR 

x{R + R 2 - z°) 



(C.8) 



,cos 2 9~ 
sin 9~ 



-n 



, sin V 

cos 6~ 



+ 1T 
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